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1. Introduction 

The study of the statistics of the fluctuating gravitional force in infinite homogeneous 
systems was pionered by Chandrasekhar & von Neumann in two classical papers (Chan- 
drasekhar & von Neumann 1942, 1943 hereafter CN43) and in several other papers by 
Chandrasekhar (1941, 1943a,b,c,d,e, 1944a and 1944b). The analysis of the fluctuating 
gravitional field, developed by the quoted authors, was formulated by means of a statis- 
tical treatment. In their papers Chandrasekhar & von Neumann considered a system in 
which the stars are distributed according to a uniform probability density, no correlation 
among the positions of the stars is present and where the number of stars constituting 
the system tends to infinity while keeping the density constant. The force F, for unit 
mass, acting on a star of a generical star system is given by the well known equation: 



where M is the mass of the i-th field star, N is the total number of the stars in the 
system and is the position vector of the i-th test star relative to the field one. The 
summation includes all the neighboring stars. The motion of the stars in the neighborhood 
of the test star produces a time variation of F. The exact dependence of F on the position 
and time cannot be exactly predicted, while it is possible to study the fluctuation of F 
from a statistical point of view. 

Two distributions are fundamental for the description of the fluctuating gravitional field: 

— 1) VK(F) which gives the probability that a test star is subject to a force F in the 
range F, F+ dF; 

— 2) W(F, f) which gives the joint probability that the star experiences a force F and 
a rate of change f, where f = dF/dt. 

The first distribution, known as Holtsmark's law (Holtsmark 1919), in the case of homoge- 
neous distribution of the stars, gives information only on the number of stars experiencing 
a given force but it does not describe some fundamental features of the fluctuations in the 
gravitional field such as the speed of the fluctuations and the dynamical friction. These 
features can be described using the second distribution W(¥, f). As shown by CN43, the 
speed of the fluctuations can be adequately expressed in terms of the mean life of a state 
F: 




(1) 



T = 



F| 



(2) 
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where (|f 2 |) is the second moment of f. Hence, for the definition of the speed of 
fluctuations and of the dynamical friction one must determine the distribution W(F, f). 
For a test star moving whith velocity v in a sea of field stars characterized by a random 
probability distribution of the velocities, $(u), we may write: 

(V) = (u) - v = -v (3) 

where V represents the velocity of a typical field star relative to the one under con- 
sideration, u denotes the velocity of a field star. This asymmetry of the distribution of 
the relative velocities produces, as shown by CN43, a deceleration of the test star in the 
direction of motion. This effect is known, after Chandrasekhar papers, as "dynamical 
friction" . Some information on dynamical friction can be obtained by means of the first 
moment of f. As shown by CN43: 



3F-v „ 
v- -r=^- F 



(4) 



where m is the mass of a field star, n; is the local density, and B({3) is definited in 
CN43 Eq. (98). These equations show that the amount of acceleration in the direction of 
—v, when v • F < 0, is greater than that in the direction +v, when v • F > 0. The star 
suffers a deceleration being the a priori probability that v • F > equal to the probability 
that v • F < 0. 

Several authors have stressed the importance of stochastic forces and in particular dynam- 
ical friction in determining the observed properties of clusters of galaxies (White 1976; 
Kashlinsky 1986, 1987) while others studied the role of dynamical friction in the orbit 
decaying of a satellite moving around a galaxy or in the merging scenario (Bontekoe & 
van Albada 1987; Seguin & Dupraz 1996; Dominguez-Tenreiro & Gomez-Flechoso 1998) 
which is not only the framework for galaxy formation picture in hierarchical cosmological 
models, but also important for the study of particular aspects of the evolution of a num- 
ber of astronomical systems, such as galactic nuclei, cD galaxies in rich galaxy clusters. 
Finally, here, we want to remember as other the statistical description of dynamical fric- 
tion other works have been made on it based on different approaches for example: Fokkcr 
Planck equation based polarization cloud (Rosenblunth et al. 1957; Binncy & Trcmaine 
1987); resonant particle interactions (Tremaine & Weinberg 1984; Weinberg 1986); fluc- 
tation dissipation (Bcrkcnstein & Maoz 1992; Maoz 1993). 

Chandrasekhar's theory (and in particular his classical formula - see Chandrasekhar 
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1943b) is widely employed to quantify dynamical friction in a variety of situations, even 
if all the theory developed by the quoted authors is based on the hypothesis that the 
stars are distributed uniformly and it is well known that in stellar systems, the stars are 
not uniformly distributed, (Elson ct al. 1987; Wybo & Dejonghe 1996; Zwart ct al. 1997) 
as well as in galactic systems, the galaxies are not uniformly distributed (Peebles 1980; 
Bahcall & Soneira 1983; Sarazin 1988; Liddle, & Lyth 1993; White et al. 1993; Strauss 
& Willick 1995). It is evident that an analysis of dynamical friction taking account of 
the inhomogeneity of astronomical systems can provide a more realistic representation 
of the evolution of these systems itself. Moreover from a pure theorethical ground we 
expect that inhomogeneity affects all the aspects of the fluctuating gravitational field 
(Antonuccio & Colafrancesco 1994; Del Popolo 1994; Del Popolo et al. 1996; Del Popolo 
& Gambera 1996, 1997; Gambera 1997). Firstly the Holtsmark distribution is no more 
correct for inhomogeneous systems. For these systems, as shown by Kandrup (1980a, b, 
1983), the Holtsmark distribution must be substituted with a generalized form of the 
Holtsmark distribution characterized by a shift of VF(F) towards larger forces when in- 
homogeneity increases. This result was already suggested by the numerical simulations of 
Ahmad & Cohen (1973, 1974). Hence when the inhomogeneity increases the probability 
that a test particle experiences a large force increases, secondly VF(F, f ) is changed by 
inhomogeneity. Consequently, the values of the mean life of a state, the first moment of 
f and the dynamical friction force are changed by inhomogeneity with respect to those 
of homogeneous systems. 

This paper must be intended as the first part of a work pointed to: 

— a) study the effects of inhomogeneity on the distribution functions of the stochastic 
forces and on dynamical friction (present paper); 

— b) test the result against N-body simulations; 

— c) find a formula that describes dynamical friction in homogeneous and inhomogeneous 
systems only on the basis of the statistical theory. 

Before continuing we want to stress that when we speak of inhomogeneity we refer to 
inhomogeneity in position distribution and not to that of velocity distribution. Our work 
follows the spirit of Kandrup's (1980a) in the sense that we are interested in the effect 
of a non-uniform distribution in position of stars on the distributions of the stochastic 
force. 

The plan of the paper is the following: in Sect. 2 we sketch the calculations needed 
to obtain the distribution function W(F,f) after having released the hypothesis of ho- 
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mogeneity. The complete calculations are developed in the appendix. Then in Sect. 3 
we calculate the first moment of f and in Sect. 4 we show how dynamical friction is 
influenced by inhomogeneity. Finally, in Sect. 5 we draw our conclusions. 

2. The distribution function W(F, f) in inhomogeneous systems 

To calculate W(F,f) in an inhomogeneous system we consider a particle moving with a 
velocity v, subject to a force, per unit mass, given by Eq. (Q) and to a rate of change 
given by 



V 4 toifaVi) 



i— 1 

where Vj is the velocity of the field particle relative to the test one. 
The expression of W(F, f) is given following Markoff's method by (CN43): 



(6) 



i pQO p CO 



'0 JO 

{exp[-i(k* + S*)]}dkdS (7) 
with A(k, S) given by 

A(k, S) = e -" c ( k - S ) (8) 
being 

/• OO /> OC f- oo 

C(k, S) = / / / r [1 - exp i(k* + £*)] drdVdM (9) 

J J ~oo J —oo 

where n is the average number of stars per unit volume while €> and are given by 
the following relations: 

M 

* = G^r (10) 



dt 



V 3r(rV) 

U3 |_|5 



(11) 



and t(V, r, M)dVdrdM is the probability that a star has velocity in the range V, V + 
dV, positions in r, r + dr and mass in M, M + dM . 
Now we suppose that r is given by: 

r=^(j 2 (M)|u| 2 ) (12) 

where a is a constant that can be obtained from the normalization condition for t, 
j a parameter (of dimensions of velocity" 1 ), ip an arbitrary function, u the velocity of 
a field star. In other words we assume, according to CN43 and Chandrasekhar & von 
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Newmann (1942), that the distribution of velocities is spherical, i.e. the distribution 
function is ip(u) = ip(j 2 (M)\u\ 2 ), but differently from the quoted papers we suppose 
that the positions are not equally likely for stars, that is the stars are inhomogeneously 
distributed in space. A lenghty calculation leads us (see Appendix for a derivation and 
the meaning of simbols) to find the function A(k, X): 

A(k,E) = e -^ V {l- igp(k,-£) 

~ -(3+p) 

+ bk^^} ■ [Q(S) + feR(E)]} (13) 

This last equation introduced into Eq. (0) solves the problem of finding the distribution 
W(F,f) and makes it possible to find the moments of f that give information regarding 
the dynamical friction. 

3. Evalutation of / 

As we stressed in the introduction, the study of the dynamical friction is possible when 
we know the first moment of f . This calculation can be done using the components of f 
(/i, fj, /fe) in the system of coordinates previously introduced. We have that: 

h ~ W(F) ( > 

and similar equations for the other components of the force. The distribution function 
W(F), giving the number of stars subject to a force F, can be calculated as follows: 

-1 /"OO />oo />0O 

w ™ = I L feMkF+Sf,,) 

•A(k, S) dkdSdf (15) 
integrating we find: 

WCF) = / { e [- a2fc(3_P)/2 ]} . ksin(kF)dk (16) 

^ Z F J 

This equation gives the generalized Holtsmark distribution obtained by Kandrup (1980a 
- Eq. 4.17) and provides the probability that a star is subject to a force F in a inhomo- 
geneous system. 

As previously stressed, to calculate the first moment of f we need only an approximated 
form for ^4(k, £): 

A(k, S) = {el -^ (3 - p)/2 ] } . [i - w(k , £)] (17) 



Using this last expression for A(k, S) and Eq. ( |A36| ), Eq. (|A37]), Eq. @, Eq. ©, 



Eq. <\VA) and performing a calculation similar to that by CN43 the first moment of f is 
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given by: 



1 \ 3_p 



A(p) ■ B{p)~ 



a~GML(/3) 



3F-v 



where 
L(/3) = 6 



i>oo 


r (3-p) " 






e W/3) * 
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COS X 



POO 


r (3-p) i 


Jo 





smi 

a;(2-p)/2 x p/2 

sinx 

ax 



2;(p-2)/2 

and for p = Eq. ( |18| ) reduces to: 
- a 



-Gm 



6 



and consequently 



[ L(J3) ' 




3F- v 




p=0 


L v |F| 2 



[i(/?)]p=0 



6 sin a; 



6 cos a; + 2x sinx 



rZ.r 



this last expression can also be written as: 

[L{fi)] P=0 = 3tt / H(0)d/3 - 7T/3H(p) 
Jo 

being 

H(/3) = —[ (e[ ( ^ )3/2 ]) -xsinxdx 
In this way we can written Eq. ( ^p| ) as: 



2tt 

Gmn 

3 



3 ■ J P g(/3)d/3 
P ■ H(P) 



1 



3Fv 



F 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



this last equation coincides with Eq. (105) by CN43. 
The results obtained by us for an inhomogeneus system are different [see Eq. ([L8[)], as 
expected, from that obtained by CN43 for a homogeneous system (CN43 - Eq. 105). At 
the same time it is very interesting to note that for p = (homogeneous system) our result 
coincides, as obvious, with the results obtained by CN43. In a inhomogeneous system, in 
a similar way to what happens in a homogeneus system, f depends on v, F and 9 (the 
angle between v and F) while differently from homogeneous systems, f is a function of 
the inhomogeneity parameter p. The dependence of f on p is not only due to the functions 
A(p), B(p) and to the density parameter a but also to the parameter P = |F|/Qjj. In fact 
in inhomogeneous systems the normal field Qh is given by Qh — GM(aB(p)/2) 2 /( 3 ~ p \ 
clearly dependent on p. 
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4. Dynamical friction in inhomogeneous systems 

The introduction of the notion of dynamical friction is due to CN43. In the stochastic 
formalism developed by CN43 the dynamical friction is discussed in terms of f : 

3F-v 



f = — ^Gmn5(/?) 
where 



(25) 



3 • fi 3 W(0)d0 , , 

Bt3 » - tmr ~ 1 (26) 

and (3 = \F\/Q H = |F|/2.603GMn 2 / 3 . As shown by CN43, the origin of dynamical 
friction is due to the asymmetry in the distribution of relative velocities. As previously 
told, if a test star moves with velocity v in a spherical distribution of field stars, namely 
</>(u) then we have that: 

V = u=v = -v (27) 
The asymmetry in the distribution of relative velocities is conserved in the final Eq. 



(|25D. In fact from Eq. (25) we have: 

m = f GMnm .f (28) 

(CN43). This means that when v • F > then > 0; while when v F < then 
— L_ — L < 0. As a consequence, when F has a positive component in the direction of v, 
|F| increases on average; while if F has a negative component in the direction of v, |F| 
decreases on average. Moreover, the star suffers a greater amount of acceleration in the 
direction — v when v • F < than in the direction +v when v • F > 0. 
In other words the test star suffers, statistically, an equal number of accelerating and 
decelerating impulses. Being the modulus of deceleration larger than that of acceleration 
the star slows down. 

At this point we may show how dynamical friction changes due to inhomogeneity. From 
Eq. ( p!s|) we see that ^ differs from that obtained in homogeneus system only for the 
presence of a dependence on the inhomogeneity parameter p. If we divide Eq. ( |l8| ) for 
Eq. (|25|) we obtain: 



,jF 



dt J 

/ Horn. 

1\& 3a^L(f3)B(p)^ -A(p) 

— — \ z v) 



n-TT 2 H{f3)B((3)(3^r 



If we consider a homogeneous system, p = 0, the previous equation reduces to: 
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dF 

dt 



Inh. 



dF 

dt 



= 1 (30) 



Horn. 

In the case of an inhomogeneous system, p ^ 0, we see that: 



dF 

dt 

Horn. 

where 



Inh. _ „p/(3-p) 



= n p/ V- p> F(0(n,p)) (31) 




Fig. 1. The function F{0) for several values of the inhomogeneity parameter p; solid line p = 0, 
dashed line p = 0.1, dotted line p = 0.5 



As we show in Fig. 1, this last equation is an increasing function of p. This means 
that for increasing values of p the star suffers an even greater amount of acceleration in 
the direction — v when v • F < than in the direction +v when v • F > 0, with respect to 
the homogeneous case. This is due to the fact that the difference between the amplitude 
of the decelerating impulses and the accelerating ones is, as in homogeneous systems, 
statistically negative, but now larger, being the scale factor greater. This finally means 
that, for a given value of n, the dynamical friction increases with increasing inhomogeneity 
in the space distribution of stars (it is interesting to note that this effect is fundamentally 
due to the inhomogeneity of the distribution of the stars and not to the density n). In 
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other words two systems having the same n will have their stars slowed down differently 
according to the value of p. This is strictly connected to the asymmetric origin of the 
dynamical friction. 

In addition, by increasing n the dynamical friction increases, just like in the homogeneous 
systems, but the increase is larger than the linear increase observed in homegeneous 
systems. 

5. Conclusion 

In this paper we calculated the distribution function and the first moment of f = dF/dt 
for an inhomogeneous system. We obtained an expression relating f = dF / dt to the degree 
of inhomogeneity in a gravitational system. In the last part of the paper we showed the 
implications of this result on the dynamical friction in a inhomogeneous system and in 
particular how inhomogeneity acts as an amplifier of the asimmetry effect giving rise to 
dynamical friction. Here we want to stress that the distribution function that we have 
obtained is valid for every inhomogeneous system and consequently it is more general 
than the distribution function obtained by CN43, that is reobtained when we assume 
p = in our model. Moreover Kandrup's (1980a) theory of the stationary distribution 
W(F) for inhomogeneous systems is reobtained in the limit t — > oo. 

Acknowledgements We are grateful to E. Recami and E. Spedicato for helpful and 
stimulating discussions during the period in which this work was performed and the 
referee Prof. P. Laguna for some useful comments and suggestions. 

Appendix: Derivation of j4(k, X) 

In order to find A(k, S) we have to calculate the integral in Eq. (||) (see Sect. 2). We 
change the integration variable in Eq. (^|) from r to $ using the following relations: 

dr = -i(GAi) L5 $- 45 # (Al) 



we have that Eq. (|^) can be rewritten as: 

-i /> DC p OC 

C(k,E) = - / / {GMf^dVdM 
2 Jo J -oo 

'GM\ 



OO 

(k*+S*) 



T 

oo 



1 -e l 



a 



j (A3) 



and if we express the product S • W as a function of $ (CN43 - Eq. 18) we can write 
C(k, S) in the following form: 
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-1 POO f oo 

G(k, E) = -G 15 / / M 15 L>(k, H)dVdM 
where D(k, £) is given by: 



D(k,£) = 

or equivalently 
D(k,E) = 



x _ e i(k*+S*) 



1 - e 



(ik€>) 



re 



;(k#) 



] • [1 - e ( jS *)] ■ $- 4 - 5 d* 



The hrst integral in Eq. (A£) can be easily calculated: 



h - 

with 
B(p) 



(ik*) 



1 -e 



z — sin z 
z (l-v)li 



47Ta (3-P) 

• A; ^ 



{GM)p/ 2 



dz 



For p = Eq. (A7) gives: 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



that coincides with the first term in the right hand side of Eq. (22) in CN43. 
The second integral is more difficult to evaluate and is given by: 

/oo 
r [ e (**)] . [i - e (* s *)] . $- 4 - 5 d* 
-OO 

/oo / jf. \ p/2 
[^[l-e^l ( Al °) 

We have obtained this last result using Eq. ([To]) and Eq. ([j~2]) (see Sect. 2). Since we are 
interested in the moments of for a given F we need only the behaviour of the function 
A(k, E) for £ => 0. To obtain this we expand the term [1 — e^^^'] in powers of £ and 
\l/ in Eq. ( A10| ) and then we have: 

(£*) 2 " 



(GM)p/2 



-j£# 



(All) 



or 



(GM) 
a 



^) p / 2 7-oo 

(S*) 2 • $(p- 9 )/ 2 . [ e C ik *)]<i* 

J — OO 



2(GAf)p/ 2 
= -Di + D 2 



(A12) 



where £ • is given by 
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£ • * = $ 1 - 5 (S 1 V) - 3$- a5 (*V) • (A13) 
being E-^ = £/(GM) 0,5 . Substituting these last equations in Eq. (A12) we have: 



£>i(k,£) 



EjV 



(GM)p/ 2 J_ 

. $ ( P -6)/2 . [ e (ik*) ]d$ 



(*V) • 



(A14) 



Now, we evalutate this integral following CN43. We introduce a system of coordinates 
with the z-axis in the direction of k and letting Si = (si, S2, S3), V = (Vj.!^2)^3) and 
1$ = (l,m,n) = (sin 9 cos w, sin 9 sin w, cos 9). The result of this integration is given by: 



Di(k,E) 
where 
A(p) = 



(GMk)p/ 2 



(31 Vi + S 2 ^2 " 2S 3 V 3 ) • A(p) 



SIM 



3 sin x 3 cos cc 
+ 



x (4-p)/2 x (S-p)/2 x (6-p)/2 



dx 



For p = Eq. ( A15 ) becomes 



£>i(k,E 



(51^1+52^2-253^3) 



(A15) 



(A16) 



(A17) 



which coincides with Eq. (37) in CN43. 



The second integral in Eq. ( A12) is: 



£> 2 (k,E 



2(GM)p/ 2 
If we define the variable 



r (s*) 2 $^ . [ e ( ik *)]d* 



(A18) 



(A19) 



and we introduce Eq. (A13) in Eq. (A18) we obtain: 

-1 



r> 2 (k,£) 



2(GM)p/ 2 



2tt 



,(izt)i 



1 JO 



[ElV - 3(1$ • V)(I« ■ E1)] 2 z^dwdtdz 



(A20) 



where i = cos# and u; is the azimuthal angle. The integral can be calculated regarding z 
and t as complex variables and using an appropriate chosen contour (see Chandrasekhar 
& von Newmann 1942). The result of the integration in terms of the original variable 
S = (E 1 ,E 2 ,E 3 ) is: 

„ ^ -™£r (3+p)/2 
r> 2 (k,53) - „ r 



2(GM)p/ 2 

IT /p 



P _ 3 
2 2 



{ C0S [- 4 (f + 1 ) 



1 



4 V3 



}{E 2 [a^ 2 + 6y 2 2 + C y 3 2 ] 



-/(S 2 S 3 y2^3 + E1E3V1V3) + eE 2 Sil/ 2 ^i} 



(A21) 



where 



13 



24 



216 



p + 1 (p+l)(p-3) (p + l)(p-3)(p-7) 
-72 



(p+l)(p-3)(p-7) 
-36 



(p-3)(p-7) 
-2 12 18 



p+1 3 — p 7 — p 

-4 48 



288 



p+1 (p+l)(p-3) (p+l)(p-3)(p-7) 
-4 24 



12 144 

+ 



p+1 (p+l)(p-3) 3-p (p-3)(p-7) 
then we have obtained: 

4-!riak~ p/2 



D(k, E) = h + h 



A 1 -S+E 

Airak 2 



B(p) - 



(GM)p/ 2 ^' (GM)p/ 2 
(EiVi + E 2 V 2 - 2£ 3 F 3 ) • A(p) + Z? 2 



Now substituting Eq. (A28) in Eq. (A4) we have 



C(k,E) = 27ra(GM) (3 - p)/2 fc( 3 - p )/ 2 B(p) - 27rm/fc- p/2 
•(GM) (2 - p)/2 (EiVi + E 2 l/ 2 - 2S3V3) • A(p) 

_™( G M)(l-p)/2 fc -(3+p)/2 r ^| + 3 
• I COS 

■{S 2 [aF^ + 6F 2 + cF 3 ] + + aF 2 + cF 3 



" 7T 








" 7T 








(I- 


— sin 




(f- 


~4 








~4 ' 







+E 2 [cF 1 + eV 2 + dF 3 ] + /(S 2 S 3 ^ 2 y 3 



+EiS 3 y 1 y 3 )+eS 2 Eil/ 2 Fi]} 



(A22) 
(A23) 
(A24) 
(A25) 
(A26) 
(A27) 



(A28) 



(A29) 



where we have used bars to indicate that the corresponding quantities have been averaged 
with the weight function r(V, M). 

If the distribution of the velocities of the field stars, u, is spherical and the test star moves 
with velocity v we have: 



V = u — v = — v 

and also: 

V u = v u 

VI = l\u\ 2 + vl 



(A30) 

(A31) 
(A32) 
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V U V V = v u v v (A33) 

and if we use the system of coordinates introduced by CN43 that is V\ — |v| serry, 
v 2 = and V3 = |v|cos7 where 7 is the angle between k and v we can simplify Eq. (A29) 
and then the Eq. (||) (see Sect. 2) becomes: 



A(k,S) 



(GM) 1 ^ |v| • [E lS in7 - 2E 3 cos 7] ■ A(p) 



a ,., i-p , -(3+p) _/p 3 
8 1 F 2+2 



•{' 



" 7T 








" 7T 




")] 




(I- 


— sin 




(f- 


~4 








~4 ' 







u 



~ (a + & + c)(E 2 + E 2 ) + ^(2c + d) 
+ |v| 2 • [(asin 2 7 + ccos 2 7)Ei + (6sin 2 7+ 
+ c cos 2 7)^2+ (csin 2 7 + dcos 2 7)E 2 + 
+ /E1S3 sin7COS7]} ] 

where a = Anna. We define, now, the following constants and functions 

~a = ~[GM)^ -B(p) 

9 = |(GM)^.|v|-4(p) 



p(k,S) 



[Ei sin 7 — 2E3COS7] 



(A34) 

(A35) 
(A36) 
(A37) 



{«* [-1 (I + 



7T 

4 



Q(S) 



^3 , 



- K+6 + C )(E 2 + E 2 ) + ^(2 C + d) 



(A38) 
(A39) 



i?(E) = (a*sin 2 7 + ccos 2 7)E 2 + (frsin 2 7 + ccos 2 7)Ej 
+ (csin 2 7 + dcos 2 7)E 2 + /E1E3 sin7 COS7 



(A40) 



so we can re- write Eq. (A34) as 



A(k,E) 



~ _(3+ p ) 

bk- L *- L ]-[Q(2) + kR(2)]} 



(A41) 
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